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ABSTRACT 

Using a simple accelerating model and an assumption that 7-rays originate from p — p 
collisions for a 7r° model, we investigate 7-ray sources like V407 Cygni in symbiotic 
stars. The upper limit of their occurrence rate in the Galaxy is between 0.5 and 5 
yr~^, indicating that they may be an important source of the high-energy 7-rays. The 
maximum energies of the accelerated protons mainly distribute around 10^"'^ eV, and 
barely reach 10^^ eV. The novae occurring in D-type SSs with ONe WDs and long 
orbital periods are good candidates for 7-ray sources. Due to a short orbital period 
which results in a short acceleration duration, the nova occurring in symbiotic star 
RS Oph can not produce the 7-ray emission like that in V407 Cygni. 
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1 INTRODUCTION 

Symbiotic stars (SSs) are usually interacting binaries, com- 
posed of a cool star, a hot component and a nebula. The 
hot component is usually a white dwarf (WD). The cool 
component is either a normal red giant (RG) in S-type or a 
Mira variable surrounded by an optically thick dust shell in 
D-type. Symbiotic novae are a small subclass of thermonu- 
clear novae which occur on a WD surface fueled by mass 
accreted from an RG. The y can pro duce sup er-soft or sof t 
X-ray emission in SSs ( Mtirset et al.i ri997: Zh u et al.ll2010l ). 

V407 Cygni is a D-t ype SS cons i sting of a Mira-type 
pulsating RG. Recently, lAbdo et al.l ((20101) reported the 
Fermi Large Area Telescope detection of variable 7-ray 
emission (0.1—10 GeV) from the nova of SS V407 Cygni. 
They explained that the 7-ray spectrum originate from 
proton-proton (p — p) interac tion by 7r° modefj], but in - 
verse Compton scattering (See lBlumenthal k, Gouldl (|l970l ) 
for an exclusive review) of infrared photons from the RG by 
electrons cannot be ruled out. The above two mechanisms 
producing 7-rays need the high-energy protons or electrons. 
Cosmic rays with high energies are thought to originate from 
supernovae remnants (SNRs). In order to having an effi- 



* E-mail: guolianglv@gmail.com (LGL) 

^ tt'^ model is that secondary neutral pion decays 7-rays from 
proton-proton (p-p) collisions, i.e., the reaction considered is 
p -\- p — > tt'^ + X and the decay is -k^ — > 27, where p repre- 
sents a proton, tt'' represents a neutral pion, and X represents 
any co mbination of particles (see the references of iKamae et al.l 

iooi)) 



cient acceleration mechanism the theory of diffusive shock 
acceleration has b een d eveloped (see a recent review from 
iMalkov fc Drurvl (|200ll )). Diffusive shock acceleration ap- 
plies only for particles with a Larmor radius larger than the 
typical shock thickness. If electrons and protons are equilib- 
rium in the shock, the Larmor radius of electron is a factor 
(me/mp)^''^ smaller than that of proton, where ttic and rrip 
are masses of an electron and a proton, respectively. Only 
electrons which are already relativistic can cross the shock 
and start acceleration. Therefore, protons are accelerated 
more easily than electrons under the mechanism diffusive 
shock acceleration. 

In this Letter, we investigate the likely possibility of the 
symbiotic novae producing 7-rays, then assume that 7-rays 
originate from p~p collisions for a vr" model, and investigate 
the 7-ray sources in SSs. 



2 ACCELERATION MODEL 

In general, the theory of diffusive shock acceleration is used 
for SNRs to explain the emission of cosmic rays with high 
energy. In this work the acceleration model used for the sym- 
biotic novae like V407 Cygni is similar with that used for 
SNRs. There is usually a low particle-density circumstellar 
medium (GSM) around SNRs. However, a symbiotic nova 
like V407 Cygni is usually embedded in a dense GSM which 
is mainly formed from stellar winds lost by the red giant. 
This difference makes it possible to efficiently accelerate pro- 
tons in the symbiotic novae. 

According to the theory of diffusive shock acceleration 
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for SNRs, the maximum attainable energy for cosmic rays 
is determined by the size of accelerator, the magnetic field 
of the CSM and the energy losses resulting from adiabatic 
processes and synchrotron processes. The siz e depends o n 
the explosion evolution of SNR. According to iKirkI (|l994l ). 
the explosion includes three phases: free-expansion phase, 
Sedov- Taylor phase and snow-plough phase. During the free- 
expansion phase, the kinetic energy of ejecta remains un- 
tapped, and the particle acceleration is not significant. Once 
the mass swept-up by the shock becomes comparable to the 
mass of the ejecta Mejc, the explosion enters the Sedov- 
Taylor phase. The acceleration efficiency is the highest in 
this phase. The maximum energy for protons is approxi- 
mately given by jSchure et al.ll201ol ') 



■iZeBVlta-v _ ZZeBV.^RsT 



(1) 



where Z is the charge number, e is the elementary electric 
charge, c is the speed of light and ^o- is a relation between the 
compression ratio of the density and magnetic field. Here, we 
assume that the magnetic field is parallel to the shock nor- 
mal, which means = 20. B is the magnetic field strength, 
and tsT ~ -Rsi/V^h is the duration for which the particles 
stay in the Sedov- Taylor phase, where Vsh is the shock speed 
and _RsT is the radius where the shock sweeps up the matter 
whose mass is equal to that of the ejecta. 

In general, for SNRs, the magnetic field B of the CSM 
is ~ /iGauss, and Rst is ^ pc due to the high mass of 
the ejecta A^a) and the low particle density of CSM 
l|McKee fc Truelovel Il995l ) . Typical tsT is several hundred 
years. However, in symbiotic novae like V407 Cygni, there 
are some physical conditions which are greatly different from 
those in SNRs: 

(i) The magnetic field of the CSM is the magn etic field of 
the st ellar wind from the RG, and it is given bv lBode et al.l 
||1985[) 



(2) 



where k is Boltzmann's constant, m = 10"'^'* g is the mean 
particle mass, and Tg is the temperature of the stellar wind. 
V407 Cygni is a D-type SS in which the RG has dust shells. 
T he temperature fo r the d ust condensation zone is ~ 1000 
K iGail fc Sedlmavil (|l999l ). Here, we take Tg = lOOOK. The 
density p is given by 



Ml 



^-KV^iB? + a2 _ 2Racosl 



(3) 



where Ml is the mass-loss rate of the RG, Vw is the stellar 
wind velocity, a is the binary separation, a distance R and 
polar angle 9 is from the WD center towards t he RG. For 
simpli city, we only consider 6 = 0°. As shown bv lBode et al.l 
l|l985l ). B is ~ 10"^ Gauss, which is 10^ times higher than 
that in the CSM of SNR. 

(ii) For a typical nova, the mass of ejecta is ~ IO^^Mq 



l Yaron et al]l2005l ). which is much less than that of the 
ejecta in a typical supernova. Furthermore, we note that 
the duration of matter ejecting in a typical nova is ~ several 
days or tens of days, and the matter ejected has an average 
expansion velocity Vkv over the whole ejecting matter phase 
and a maximal expansion velocity Knax. This means that a 
part of the matter ejected has the high expansion velocity 
Vinax- The shock in a nova is mainly produced by the matter 



ejected with high velocity. We assume that Vsi-i = Vmax, and 
use a parameter 77 to define a ratio of the mass ejected with 
a velocity of Knax to the whole ejecta. Rst can be given by 
the following equation: 

AixR^pAR (4) 

-'WD 

where 7?wd is the radius of WD. lAbdo et al.l (|2010l ) found 
that the peak fiux in 7-rays was observed after 3-4 days of 
a nova outburst from V407 Cygni on 10 March 2010. This 
implies that tsT in V407 Cygni should shorter than 3 days. 
In our model tsT depends on the parameter r). We find that 
tsT ~ days when -q ^ 0.01. 

The novae occurring in SSs are surrounded by the dense 
stellar winds from the RGs. They offers an environment for 
high efficient particle acceleration. Therefore, they may be 
an important source of the high-energy 7-rays in the Galaxy. 



3 SYMBIOTIC STARS 

In general, SSs are the detached interacting binaries in which 
the WDs accrete the matter of the RGs via stella r winds. 
By a population synthesis method, iLii et all (| 20061 ) carried 
out a detailed investigation of SSs. They found that the 
occurrence rate of the novae in SSs is greatly affected by 
common-envelope ev olution and the stella r wind velocity K v 
of the RG. Following iLii et all (|2006l ') and lZhu et al] |20lO), 
for common-envelope evolution in different simulations we 
use an acoAco = 0.5 in a-algorithm and 7 = 1.75 in a 7- 
algorithm, respectively; for the stellar wind, = I'^csc 
where v^sc is the escape velocity and Vw is determined by 
the relation between th e mass-loss rates and the terminal 
wind velocities fitted by I Winters et all (|2003l ) as: 

logio(M/MQyr-i) = -7.40 + ^ logio(Vw/kms-i). (5) 

In this work we consider three cases with different input 
parameters: 

(i) in case 1, QccAcc = 0.5 and 1/„ = \vcsc\ 

(ii) in case 2, 7 = 1.75 and Vw ~ \vcsc\ 

(iii) in case 3, aceAco = 0.5 and Vw taken as Eq. ([5]). 
Using t he m odel of SSs and grid for novae in 

lYaron et al.l (|2005l ). we can estimate -Emax in which 6 = Q° 
and 17 = 0.01 for every nova. According to iKamae et al.l 
(|2006l ) , p — p interaction can occur when the energy of pro- 
ton is higher than 10® eV, which results in 7-ray emission. 
Therefore, we assume that the novae in SSs are 7-ray sources 
ff the £P ax in Eq. (O is higher than 10® eV. 



4 RESULTS 



Using a pop ulation synthesis method described in lLii et al.l 
(|2006l . bOOSl ). we model 10® binary systems which gives a 
statistical error for our Monte Carlo simulation lower than 
5 percent for the symbiotic novae. In order to estimate the 
occurrence rate of the 7-ray sources like V407 Cygni, we 
assume one binary with primary mass more massive than 0.8 
Mq is formed annually in the Galaxy. We do not consider 
the energy losses of the accelerated protons via adiabatic 
process and synchrotron process. Therefore, we overestimate 
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Figure 1. Gray-scale maps of WD's masses vs. orbital periods 
for SSs as 7-ray sources in cases 1, 2 and 3. The gradations of 
gray-scale correspond to the regions where the number density of 
systems is, respectively, within 1 - 1/2, 1/2 - 1/4, 1/4 - 1/8, 1/8 
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Figure 2. Similar to Figure[T] but for the maximum energy 
of the protons accelerated vs. the shock velocity Vgij. 



the maximum energy of the proton and the occurrence rate 
of the 7-ray sources in this work. 

We select symbiotic novae as 7-ray sources if iS^ax of 
the accelerated protons is larger than 10^ eV. Our model 
shows that the upper limits of the occurrence rates of 7-ray 
sources like V407 Cygni in SSs are 0.5 yr~^ in case 1, 2.0 
yr~^ in case 2 and 5.0 yr"*^ in c a se 3, respectively. Com- 
pared with the results in iLii et all (|2006l ). about 15% of the 
novae in SSs for cases 1 and 2 can produce 7-ray emission, 
and it is 40% in case 3 because the wind velocity in Eq.([5ll is 
favorable for a strong nuclear outburst. If the Galactic cos- 
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Figure 3. Similar to Figure [T] but for the duration of Sedov- 
Taylor phase fsT vs. the magnetic field of the stellar wind from 
RGs. 



mic rays originate from the supernova whose occurrence in 
the Galaxy is ~ 0.01 yr~^, we suggest that symbiotic novae 
like V407 Cygni may be another important source of the 
high-energy 7-rays. However, the contribution of the sym- 
biotic novae to total cosmic rays cannot be known until the 
spectra-energy distribution of the 7-rays from the novae is 
calculated, which will be carried out in further work. 

Figure [1] shows the distribution of the WD's masses vs. 
orbital periods for SSs as 7-ray sources. Majority of WDs in 
these SSs a re ONe WDs an d they have masses larger than 
1.3 Mq. As lLii et al.l (|2008l ) mentioned, the most significant 
property of novae occurring on the surface of ONe WDs is 
high neon abundance in the ejected materials. V407 Cygni 
may offer a chance to investigate the thermonuclear out- 
bursts on the surface of ONe WDs. However, to our knowl- 
edge, there is no observational data on the neon abundance 
of the ejecta in this nova from V407 Cygni. The peak of 
orbital-period distribution is around ~ 10* d ays in cases 1 
and 3, and it is around ~ 10'^ days in case 2. iMunari et al.l 
( .1990. ) suggested that the V407 Cygni has an orbital period 
of 43 years, which is consistent with our results. Novae in SSs 
with long orbital periods mean that RGs have high mass- 
loss rates. In our work, the majority of th e RGs have mass- 
loss ra tes higher than 5 x 1O~^M0 vr~^. [Fbrrarotti fc Gaij 
(|2006l) suggested that RGs produce significant dust when 
their mass- loss rates are higher than 3 x IO'^^Mq yr^^. 
Therefore, we consider that most of SSs as 7-ray sources 
have massive WDs, long orbital periods and RGs with dust 
shells. This means that the novae occurring in D-type SSs 
with ONe WDs are good candidates for 7-ray sources. 

RS Oph is a symbiotic recurrent nova which had pre- 
viously undergone recorded outbursts in 1898, 1933, 1958, 
1967 and 1985. It comprises a red giant star in a 455.72±0.83 
day orbital period with a white dwarf (WD) of mass near 
the Ghandrasekhar limit. Recently, RS Oph was observed 
to be undergoin g an outburst On 2006 Februry 12.83 UT 
(|Hirosowall2006 lV During the first 3 days a hard X-ray emis- 
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sion (14-25 keV) was clearly detected, and there was a 
weak detection in the 25-50 keV ba nd immediately follow- 
ing the outburst (jBode et al.ll200^ . After that, the X-ray 
spectrum in RS Oph was seen to evolve from the relatively 
h ard to a super-soft s ource state (Evolutionary details are 
in lNelson et al.l l)2008h l. However, there is no conclusive ev- 
idence for gamma-ray emission like V407 Cygni in the 2006 
outburst. The orbital period of V407 Cygni is ~ 40 times 
that of RS Oph, which means that the density of CSM 
around the WD in the later is higher than 100 times of 
that in the former if is comparable in the two binariefl 
Therefore, the duration of the diffusive shock acceleration 
in 2006 outburst of RS Oph, tsT, is too short so that pro- 
tons and electrons can not be accelerated enough energy to 
produce 7-ray emission. 

Figure[2]gives the distribution of -Bmax vs. V^h. The peak 
of £^max is at ~ 10" eV, and B^ax hardly reaches 3 x lO^'^ 
eV which is the knee of the cosmic ray spectra. If the 7- 
ray energy originating from p — p interaction is comparable 
to -Emaxj they should be in the low- frequency part of high- 
energy cosmic rays. These nuclear outbursts are very strong 
so that the ejecta have high velocity. As Figure [2] shows, Vsh 
in the symbiotic n ovae is ~ several 10 ^ km/s. According to 
the calculations of lYaron et al.l (|2005l l. in these strong nu- 
clear outbursts most of the accreted mater is expelled and in 
some cases even an erosion of the WD occurs. Therefore, the 
massive WDs in 7-ray sources do not explode as supernovae. 

Figure[3]shows the distribution of the duration of Sedov- 
Taylor phase tsT vs. the magnetic field of the stellar wind 
from RGs. The peaks of the magnetic field distribution are 
around ^ 10~^ Gauss in cases 1 and 2, and it is around 
~ 10~^ Gauss in case 3 because a high mass-loss rate results 
in high Vw (See Eq. (O). The magnetic fields of the stellar 
winds around the novae in SSs are 10^ times higher than 
those of CSM around SNRs. The peak of tsT is at ~ 3 days 
in cases 1 and 2, while it is around ~ 10 days in case 3 due 
to a high which results in the low density of stellar wind. 
However, as mentioned in the tsT are greatly affected by 
an uncertain parameter rj. 



5 CONCLUSIONS 

In this Letter, we use a toy model to investigate the 7-ray 
sources which originate from p — p collisions in the novae 
of SSs. The symbiotic novae occurring on the surface of ac- 
creting WDs are surrounded by the dense stellar winds from 
the RGs. They offers an environment for highly efficient par- 
ticle acceleration. We estimate that the upper limit of the 
occurrence rate of 7-ray sources in SSs in the Galaxy is be- 
tween 0.5 and 5 yr~^. Therefore, they may be an important 
source of the high-energy 7-rays. The maximum energies of 
the accelerated protons mainly distribute around 10^^ eV, 
and barely reaches 10^^ eV. If the 7-ray energy originating 
from p — p interaction is comparable to i5£iax, they should 



be in the low-frequency part of high-energy cosmic rays. In 
SSs as 7-ray sources, majority of WDs are ONe WDs and 
have masses larger than I.3M0, and most of RGs have dust 
shells. The novae occurring in D-type SSs with ONe WDs 
are good candidates as 7-ray sources. Due to the short or- 
bital period, the nova occurring in RS Oph hardly produces 
7-ray emission like that in V407 Cygni. 
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iMunari et al. estimated that the mass-loss rate of the 

RG In V407 Cygni Is ~ 6 x 10"'^ M© yr~^. Considering no thick 
dust shells around RS Oph, we assume that the mass-loss rate of 
the RG In R S Oph Is ~ W~'^ Mp, yr ~^. However, due to the dust- 
driven wind llGall fc Sedlmavi1ll986h , the stellar wind velocity Kv 
In V407 Cygni Is higher than that In RS Oph 



